Super-lattice, rhombus, square, and hexagonal standing waves 
in magnetically driven ferrofluid surface 
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Standing wave patterns that arise on the surface of ferrofluids by 
(single frequency) parametric forcing with an ac magnetic field are 
investigated experimentally. Depending on the frequency and ampli- 
tude of the forcing, the system exhibits various patterns including a 
superlattice and subharmonic rhombuses as well as conventional har- 
monic hexagons and subharmonic squares. The superlattice arises 
in a bicritical situation where harmonic and subharmonic modes col- 
lide. The rhombic pattern arises due to the non-monotonic dispersion 
relation of a ferrofluid. 
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The phenomena of self-organization in spatially extended 
nonequilibrium systems have attracted much attention lately 
in the scientific community world-wide [ffj. One of the re- 
cent scientific issues in this regard is to investigate the system 
in which a finite number of unstable (spatial and temporal) 
modes interact each other, while the past studies have mostly 
focussed on the limiting cases of either pattern forming sys- 
tems where a single mode is excited or systems that are highly 
turbulent with many modes. Although the mode-coupling in- 
teractions have been a long standing issue in various fields 
ranging from material science to plasma physics, their signifi- 
cance in nonequilibrium pattern formation is only begun to be 
understood. 

In a recent experiment of parametrically driven surface 
waves, Edwards and Fauve used a two-frequency excitation 
scheme and were able to generate a quasi-crystalline standing 
wave pattern with 12-fold symmetry [0]. In subsequent ex- 
periments on a similar experimental setup, Kudrolli et al. ob- 
served two different types of superlattice pattern ("type-I" 
and "type-II") in addition to the 12-fold symmetric quasi- 
crystalline pattern More recently, Arbell and Fineberg 
report yet another kinds of superlattice pattern ("SSS" and 
"2MS") [fil]. All these unusual patterns are obtained in two- 
frequency forced Faraday experiments in which two exter- 
nally excited modes interact (i. e. near bicriticality). Very 
recently, however, Wagner et al. have demonstrated that such 
bicritical condition can be achieved even by a single frequency 
forcing for a viscoelastic liquid |Q| or for a normal Newto- 
nian fluid in an extreme driving condition (very shallow filling 
depth and large shaking elevations) [||. New types of super- 
lattice pattern are revealed in those systems. The existence 
of these unusual patterns all together suggests a large class of 
previously unanticipated spatio-temporal patterns with multi- 
ply interacting modes. 






FIG. 1. Regular patterns observed on the surface of a driven fer- 
rofluid in a cylindrical container: (a) harmonic hexagon (/ = 6.0 
Hz), (b) subharmonic superlattice (/ = 6.7 Hz), (c) subharmonic 
square (/ = 11.0 Hz), and (d) subharmonic rhombus (/ = 20.0 
Hz). The total magnetic field is H (t) = Hq + AH sin(27r/t), where 
H /H c and AH/H C are fixed at 0.94 and 0.22, respectively. The 
diameter of the viewing area is 60 mm. 



In the present paper, we take a step in this direction by 
experimental study on a vertically oscillating ferrofluid that 
is driven by a sinusoidal magnetic field. There are two rea- 
sons for choosing the new system over previously investi- 
gated ones. First of all, the driven ferrofluid system can 
be easily brought near to a bicritical situation with a single 
frequency forcing - the harmonic mode as well as subhar- 
monic mode can be generated by an ac magnetic field with 
a small driving frequency (< 30 Hz) and a small magnetic 
field (< 100 gauss) [Q]. This contrasts with the case of me- 
chanically driven Newtonian fluids in which either a complex 
two-frequency forcing scheme or an extreme driving condi- 
tion is necessary to achieve a bicritical situation [g]. The other 
reason is that ferrofluid exposed to a magnetic field can ex- 
hibit a non-monotonic dispersion relation [gp. Due to the non- 
monotonicity the neural stability curve for the excitation of 
standing waves can have multiple minima JlO| , pl] ], thus, sev- 
eral spatial modes can be excited simultaneously. 

As expected, the driven ferrofluid system is quite rich to re- 
veal a variety of standing wave patterns as shown in Fig. [l| 
Two unusual patterns, (1) superlattice with a mixed-mode 
oscillation and (2) subharmonic rhombic pattern, as well as 
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usual subharmonic square and harmonic hexagonal patterns 
are observed. The observed superlattice is similar but differ- 
ent from any of the previously observed ones in its shape and 
temporal evolution. The rhombic pattern is the first example 
of its kind ever demonstrated in experiments on nonequilib- 
rium pattern formation. 

The ferrofiuid is a colloidal suspension of magnetic pow- 
der stabilized by screened electrostatic repulsion [^j. We use 
a commercially available ferrofiuid (EMG901, Ferrofiuidics) 
[density p = 1.53 gm/ml, surface tension a = 29.5 g/s 2 , 
initial magnetic susceptibility \ = 3.00, magnetic saturation 
M s = 600 gauss, dynamic viscosity r\ = 10 cp, yielding a 
critical field of the static Rosensweig-instability H c = 79.73 
gauss]. A cylindrical Teflon container containing the fer- 
rofiuid (depth = 5 mm, diameter = 85 mm) is placed in the 
center of a pair of Helmholtz-coils (Cenco instruments) with 
a inner (outer) diameter of 100 mm (300 mm). The magnetic 
field is monitored by means of a hall probe (F. W. Bell Inc., 
Model 4048), and the spatial variation of the field strength 
from the center to the outer rim is within 3%. An ac signal 
is generated from a home-built synthesizer-board and ampli- 
fied by a linear amplifier driving the total magnetic field of 
H{t) = H + AHsm(2irft). H Q is the static field, AH is 
the amplitude of ac component, and / is the driving frequency. 
Hq is fixed at 0.94ff c for all cases, while AH and / are used 
as a control parameter. 

The fluid surface is illuminated by two arrays of concen- 
tric LED rings (diameter 100 mm and 120 mm, respectively) 
located 230 mm above the surface. The flat surfaces either 
above or below the level of the surrounding fluid appear white, 
while the non-flat surfaces that scatter the light away from the 
camera appear black. The patterns are imaged using a high 
speed 256 x 256 pixel CCD camera (Dalsa, CA-D6-0256W) 
located 370 mm above the surface with a fast frame grabber 
(Matrox, Meteor-II). Maximum frame acquisition rate is 955 
frame/s, enabling us to photograph distinct phases of moving 
pattern. 

The phase diagram shown in Fig. || summarizes our results. 
Upon increasing sequence of AH, the flat surface can be un- 
stable to various crystalline states depending on the value of 
driving frequency /. Squares and rhombuses occur over a 
wide range of parameter space, whereas hexagons and super- 
lattice form only in a narrow range at low driving frequency. 
Spatio-temporal chaotic state (STC) gradually forms with in- 
creasing number of defects as AH increases. Essentially, the 
same phase diagram is obtained with a square container. 

We first focus on the sequence of bifurcations that occur 
with increasing value of / for a fixed value of AH/H C = 
0.22. The flat fluid surface becomes unstable first to a har- 
monic mode for a small /. Since the up-down symmetry is 
broken for a harmonic mode, quadratic interaction is the most 
significant nonlinearity [jlj], producing hexagonal standing 
wave pattern that is synchronous to the ac driving [Fig. |l](a) 
and Fig. |3^a)]. With a gradual increase in /, the wave num- 
ber (fcf/) increases slightly [see Fig. ^|(e)] but the hexagonal 
symmetry does not change until / = 6.7 Hz at which the sys- 
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FIG. 2. Phase diagram of standing wave patterns revealed in a 
magnetically driven ferrofiuid with Hq = 0.94H C : flat surface (FL), 
hexagon (H), superlattice (SU), square (SQ), and rhombus (RH). The 
thin shaded area between SQ and RH is bistable (i. e. SQ and RH 
coexist). Spatio-temporal chaotic state (STC) is generally observed 
at large AH. Crossing the dashed line toward the STC state, the 
crystalline states start to acquire defects. 



(a) ■ 

1 k H 
1 


(b) 

■ 

iv S ub ■ 

■ N ■ 

\ 

\ 

■ 1 


(c) 

k SQ 

■ ■ 
_ 


KRH2 

|^J<RH1 

■ 



(e) 220 



180 



[1/m] 




140 



100 



60 



6.20 6.60 7.00 7.40 



^1 



»RH1 



»RH2 



sub 



4.0 



8.0 



16.0 



20.0 



12.0 

f[Hz] 

FIG. 3. Fourier transforms of the regular patterns shown in Fig. [lj 
(a) hexagon, (b) superlattice, (c) square, and (d) rhombus. The mag- 
nitude of wavevectors associated with different patterns vs f are 
shown in (e). The insets shows the details of the transition among 
hexagon, superlattice, and square. 
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FIG. 4. Time evolution of subharmonic patterns during one driv- 
ing cycle: (a) superlattice, (b) square, and (c) rhombus. The patterns 
after one driving period (r) are shown in the bottom row. The dashed 
cross lines are placed to indicated the same location. 

tern enters into the bicritical region. At the onset of bicritical- 
ity, 6 subharmonic peaks (fc su b = ku/ v3) suddenly appear 
forming a superlattice pattern together with the primary 6 har- 
monic peaks [see Fig. |l](b) and Fig. |3|(b)]. The transition from 
hexagon to superlattice is very sharp and appears to be first or- 
der, but the exact nature of the transition could not be resolved 
within the accuracy of our current experimental setup. 

The observed superlattice pattern is similar but different 
from the ones previously observed by others [B-§]. In par- 
ticular, the snapshot of superlattice [Fig. 0(a)] taken at t = 
is basically tiled by three different cells differing in sizes and 
shapes. During one driving cycle the cell "M" returns to the 
same "M", while "S" and "L" swap each other. Upon com- 
pleting the subsequent driving cycle, "M" returns to the same, 
but "S" and "L" swap again. In other words, the superlattice 
in Fig. |l](b) has a mixed mode oscillation. This is in marked 
contrast with the purely subharmonic hexagonal superlattices 
observed by Wagner et al. that alternates between two comple- 
mentary hexagonal lattices composed with cells of same type 
[H], and also with the "Type-II superlattice" observed by Ku- 
drolli et al. that undergoes a complex stripe modulation during 
the course of oscillation [||]. 

Since the superlattice forms in a narrow region between har- 
monic hexagon and subharmonic square, it is natural to be- 
lieve that the superlattice forms by an interplay between the 
two different modes. Since both harmonic and subharmonic 
modes are present, dominant nonlinearity of inter-mode in- 
teraction is quadratic generating new subharmonic modes. 
In addition, translational invariance requires that the newly 
generated wavevectors are linear combinations of kn and 
ksQ- They would in turn interact with the dominant harmonic 
hexagonal basis, and only those which satisfy the (spatial 
and temporal) resonance conditions would survive [ |l2[ ] . Two 
out of four subharmonic square modes ksq and four out of 
new subharmonic modes satisfy the resonance conditions, and 
form a hexagonal sublattice with magnitude fc su b = kji/ Vo, 
thus contributing to the superlattice. The present mechanism 
for the generation of the superlattice is very similar to that of 
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As / is increased further, the superlattice undergoes a hys- 
teretic transition to a subharmonic square lattice; upon in- 
creasing (decreasing) sequence of /, the superlattice (square) 
transforms to the square (superlattice) at / = 7.10 Hz (/ = 
7.20 Hz). At the onset of transition to square lattice, ksQ is 
approximately kjj / v2- In the bistable regime, mixed states of 
square domains and superlattice domains are often seen. The 
observed square patterns are subharmonic just like the ones 
that are popularly seen in conventional Faraday experiments: 
the pattern alternates between two dual square lattices at every 
driving period [see Fig. |Jb)]. The ksq increases significantly 
as / increases. 

Upon increasing / further, the square lattice undergoes a 
hysteretic transition to a stable rhombic pattern with a well 
defined 2-fold symmetry [see Fig. ||(d) and (e)]. At the tran- 
sition, the ksQ branch continues to kani (primary rhombic 
basis vector) and a new basis vector knH2 that is the differ- 
ence between two neighboring kuni^ is spontaneously cre- 
ated. The oblique angle between kani and kRH2 character- 
izes a particular rhombic pattern. For the example shown in 
Fig. |3](d), knuilkRHi — 0.696 and the angle between the 
two neighboring basis wavevectors is 69.7°. The oblique an- 
gle changes but slightly as the driving frequency changes (less 
than 2% for the frequency range 17.5 - 21.0 Hz). The system 
is bistable for the range (17.5Hz</<18.0 Hz). Except for 
the bistable region, the rhombic lattice is globally stable. Fig- 
ure [|(c) clearly demonstrate that the observed rhombic pattern 
is temporally subharmonic just like the square pattern shown 
m Fig.gb). As far as we know, the observed rhombic pattern 
is the first laboratory demonstration of its kind ever revealed 
& 

At this point, we note that a recent theoretical study has 
shown that parametrically driven surface waves on ferroflu- 
ids can be unstable to multiple spatial modes simultaneously 
Jic|]. Here, the existence of multiple unstable modes origi- 
nates from the non-monotonic dispersion relation of ferroflu- 
ids. The study has shown that in one-dimensional geometry 
one can realize a situation that different domains with differ- 
ent wavevectors coexist. Such phenomenon is also verified re- 
cently in an experiment conducted in a quasi one-dimensional 
annular geometry [O]. The existence of rhombic pattern in 
our system, however, suggests that in two-dimensional space 
the (two) unstable modes rather form a deformed crystalline 
structure (i. e. rhombic lattice), instead of domains with differ- 
ent wavevectors. In our experiments, coexisting domains with 
different wavevectors are only observed in bistable regions. In 
order to check the feasibility of this hypothesis, we have ex- 
tended the earlier theoretical study by Raitt and Riecke [ |To| ] to 
two-dimensional space and have found a parameter regime in 
which well defined rhombic patterns are globally stable [M|. 
Recently, Lifshitz and Petrich have also demonstrated the ex- 
istence of a stable rhombic pattern in a generic model system 
containing two unstable spatial modes jl5|]. 

In the following, the transition to spatio-temporal chaotic 
state (STC, see Fig |^) from crystalline states is qualitatively 
described. As AH increases for a fixed value of /, the crys- 
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FIG. 5. Transitions to disordered patterns from (a) square and (b) 
rhombus, as AH increases. The frequency / is fixed at 15.0 Hz for 
(a) and 21.0 Hz for (b). AH/H C is 0.28, 0.39, 0.45 (from left to 
right) for (a), and 0.28, 0.37, 0.45 for (b). 

talline patterns gradually become disordered. Figure || shows 
two such sequences, one for the case of square and the other 
for the case of rhombic pattern. Dislocations and defects start 
to form in a square lattice near AH/H C = 0.28 [Fig. |]^a), 
left image]. As AH increases, the number of defects in- 
creases (middle image) and eventually the whole lattice be- 
comes very loose as in a typical melting process. The last 
image in Fig. ||(a) shows a state very near to STC, showing a 
very little trace of square lattice. A similar melting process is 
also observed for the case of rhombic lattice [Fig. j5|(b)]. The 
gradual deformation toward the state of STC is not limited to 
the square and rhombus but occurs in the whole range of / 
that we have explored. This observation is consistent with the 
earlier Faraday experiments conducted with Newtonian fluids 

& 

In conclusion, superlattice and rhombus, as well as square 
and hexagonal standing waves, are observed on the surface of 
ferrofluids driven parametrically with an ac magnetic field. It 
is argued that the superlattice pattern is formed by the reso- 
nance of newly generated subharmonic modes with the har- 
monic hexagonal mode already present. Due to the non- 
monotonic dispersion relation of the ferrofluids, multiple spa- 
tial modes can become unstable simultaneously with a single 
frequency forcing. As a consequence, stable rhombic patterns 
are possible and observed in the experiment. Full understand- 
ing of the observed patterns and their mutual boundaries, how- 
ever, pose quite a challenge. 

We are indebted to W. G. Choe for many helpful discus- 
sions. This work is supported by Creative Research Initiatives 
of the Korean Ministry of Science and Technology. 
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